5-Hydroxymethylcytosine in Cell-Free DNA Predicts Immunotherapy Response in Lung Cancer.

Immune checkpoint inhibitors (ICIs) drastically improve therapeutic outcomes for lung cancer, but accurately predicting individual patient responses to ICIs remains a challenge. We performed the genome-wide profiling of 5-hydroxymethylcytosine (5hmC) in 85 plasma cell-free DNA (cfDNA) samples from lung cancer patients and developed a 5hmC signature that was significantly associated with progression-free survival (PFS). We built a 5hmC predictive model to quantify the 5hmC level and validated the model in the validation, test, and control sets. Low weighted predictive scores (wp-scores) were significantly associated with a longer PFS compared to high wp-scores in the validation [median 7.6 versus 1.8 months; p = 0.0012; hazard ratio (HR) 0.12; 95% confidence interval (CI), 0.03-0.54] and test (median 14.9 versus 3.3 months; p = 0.00074; HR 0.10; 95% CI, 0.02-0.50) sets. Objective response rates in patients with a low or high wp-score were 75.0% (95% CI, 42.8-94.5%) versus 0.0% (95% CI, 0.0-60.2%) in the validation set (p = 0.019) and 80.0% (95% CI, 44.4-97.5%) versus 0.0% (95% CI, 0.0-36.9%) in the test set (p = 0.0011). The wp-scores were also significantly associated with PFS in patients receiving single-agent ICI treatment (p < 0.05). In addition, the 5hmC predictive signature demonstrated superior predictive capability to tumor programmed death-ligand 1 and specificity to ICI treatment response prediction. Moreover, we identified novel 5hmC-associated genes and signaling pathways integral to ICI treatment response in lung cancer. This study provides proof-of-concept evidence that the cfDNA 5hmC signature is a robust biomarker for predicting ICI treatment response in lung cancer.

Iterative oxidation by TET1 is required for reprogramming of imprinting control regions and patterning of mouse sperm hypomethylated regions.

Ten-eleven translocation (TET) enzymes iteratively oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxylcytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during mammalian germline reprogramming remains unresolved due to the inability to decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 (Tet1-HxD) and TET1 that stalls oxidation at 5hmC (Tet1-V). Tet1 knockout and catalytic mutant primordial germ cells (PGCs) fail to erase methylation at select imprinting control regions and promoters of meiosis-associated genes, validating the requirement for the iterative oxidation of 5mC for complete germline reprogramming. TET1V and TET1HxD rescue most hypermethylation of Tet1-/- sperm, suggesting the role of TET1 beyond its oxidative capability. We additionally identify a broader class of hypermethylated regions in Tet1 mutant mouse sperm that depend on TET oxidation for reprogramming. Our study demonstrates the link between TET1-mediated germline reprogramming and sperm methylome patterning.

Sequence based model using deep neural network and hybrid features for identification of 5-hydroxymethylcytosine modification.

RNA modifications are pivotal in the development of newly synthesized structures, showcasing a vast array of alterations across various RNA classes. Among these, 5-hydroxymethylcytosine (5HMC) stands out, playing a crucial role in gene regulation and epigenetic changes, yet its detection through conventional methods proves cumbersome and costly. To address this, we propose Deep5HMC, a robust learning model leveraging machine learning algorithms and discriminative feature extraction techniques for accurate 5HMC sample identification. Our approach integrates seven feature extraction methods and various machine learning algorithms, including Random Forest, Naive Bayes, Decision Tree, and Support Vector Machine. Through K-fold cross-validation, our model achieved a notable 84.07% accuracy rate, surpassing previous models by 7.59%, signifying its potential in early cancer and cardiovascular disease diagnosis. This study underscores the promise of Deep5HMC in offering insights for improved medical assessment and treatment protocols, marking a significant advancement in RNA modification analysis.

NSUN5/TET2-directed chromatin-associated RNA modification of 5-methylcytosine to 5-hydroxymethylcytosine governs glioma immune evasion.

Malignant glioma exhibits immune evasion characterized by highly expressing the immune checkpoint CD47. RNA 5-methylcytosine(m5C) modification plays a pivotal role in tumor pathogenesis. However, the mechanism underlying m5C-modified RNA metabolism remains unclear, as does the contribution of m5C-modified RNA to the glioma immune microenvironment. In this study, we demonstrate that the canonical 28SrRNA methyltransferase NSUN5 down-regulates ß-catenin by promoting the degradation of its mRNA, leading to enhanced phagocytosis of tumor-associated macrophages (TAMs). Specifically, the NSUN5-induced suppression of ß-catenin relies on its methyltransferase activity mediated by cysteine 359 (C359) and is not influenced by its localization in the nucleolus. Intriguingly, NSUN5 directly interacts with and deposits m5C on CTNNB1 caRNA (chromatin-associated RNA). NSUN5-induced recruitment of TET2 to chromatin is independent of its methyltransferase activity. The m5C modification on caRNA is subsequently oxidized into 5-hydroxymethylcytosine (5hmC) by TET2, which is dependent on its binding affinity for Fe2+ and α-KG. Furthermore, NSUN5 enhances the chromatin recruitment of RBFOX2 which acts as a 5hmC-specific reader to recognize and facilitate the degradation of 5hmC caRNA. Notably, hmeRIP-seq analysis reveals numerous mRNA substrates of NSUN5 that potentially undergo this mode of metabolism. In addition, NSUN5 is epigenetically suppressed by DNA methylation and is negatively correlated with IDH1-R132H mutation in glioma patients. Importantly, pharmacological blockage of DNA methylation or IDH1-R132H mutant and CD47/SIRPα signaling synergistically enhances TAM-based phagocytosis and glioma elimination in vivo. Our findings unveil a general mechanism by which NSUN5/TET2/RBFOX2 signaling regulates RNA metabolism and highlight NSUN5 targeting as a potential strategy for glioma immune therapy.

Association of tet methylcytosine dioxygenase 2 and 5-hydroxymethylcytosine in endometrioid adenocarcinoma and its clinical significance.

BACKGROUND: Aberrant DNA methylation is a vital molecular alteration commonly detected in type I endometrial cancers (EC), and tet methylcytosine dioxygenase 2 (TET2) and 5-hydroxymethylcytosine (5hmC) play significant roles in DNA demethylation. However, little is known about the function and correlation of TET2 and 5hmC co-expressed in EC. This study intended to investigate the clinical significance of TET2 and 5hmC in EC. METHODS: The levels of TET2 and 5hmC were detected in 326 endometrial tissues by immumohistochemistry, and the correlation of their level was detected by Pearson analysis. The association between the levels of TET2 and 5hmC and clinicopathologic characteristics was analyzed. Prognostic value of TET2 and 5hmC was explored by Kaplan-Meier analysis. The Cox proportional hazard regression model was used for univariate and multivariate analyses. RESULTS: Based on the analysis results, TET2 protein level was positively correlated with 5hmC level in EC tissues (r = 0.801, P < 0.001). TET2+5hmC+ (high TET2 and high 5hmC) association was significantly associated with well differentiation, myometrial invasion, negative lymph node metastasis, and tumor stage in EC. Association of TET2 and 5hmC was confirmed as a prognostic factor (HR = 2.843, 95%CI = 1.226-3.605, P = 0.007) for EC patients, and EC patients with TET2-5hmC- level had poor overall survival. CONCLUSIONS: In summary, the association of TET2 and 5hmC was downregulated in EC tissues, and may be a potential poor prognostic indicator for EC patients. Combined detection of TET2 and 5hmC may be valuable for the diagnosis and prognosis of EC.

Construction of a glycosylation-mediated fluorescent biosensor for label-free measurement of site-specific 5-hydroxymethylcytosine in cancer cells with zero background signal.

BACKGROUND: 5-hydroxymethylcytosine (5hmC) as an epigenetic modification can regulate gene expression, and its abnormal level is related with various tumor invasiveness and poor prognosis. Nevertheless, the current methods for 5hmC assay usually involve expensive instruments/antibodies, radioactive risk, high background, laborious bisulfite treatment procedures, and non-specific/long amplification time. RESULTS: We develop a glycosylation-mediated fluorescent biosensor based on helicase-dependent amplification (HDA) for label-free detection of site-specific 5hmC in cancer cells with zero background signal. The glycosylated 5hmC-DNA (5ghmC) catalyzed by ß-glucosyltransferase (ß-GT) can be cleaved by AbaSI restriction endonuclease to generate two dsDNA fragments with sticky ends. The resultant dsDNA fragments are complementary to the biotinylated probes and ligated by DNA ligases, followed by being captured by magnetic beads. After magnetic separation, the eluted ligation products act as the templates to initiate HDA reaction, generating abundant double-stranded DNA (dsDNA) products within 20 min. The dsDNA products are measured in a label-free manner with SYBR Green I as an indicator. This biosensor can measure 5hmC with a detection limit of 2.75 fM and a wide linear range from 1 × 10-14 to 1 × 10-8 M, and it can discriminate as low as 0.001% 5hmC level in complex mixture. Moreover, this biosensor can measure site-specific 5hmC in cancer cells, and distinguish tumor cells from normal cells. SIGNIFICANCE: This biosensor can achieve a zero-background signal without the need of either 5hmC specific antibody or bisulfite treatment, and it holds potential applications in biological research and disease diagnosis.

Epigenetic landscape of 5-hydroxymethylcytosine and associations with gene expression in placenta.

5-hydroxymethylcystosine (5hmC), is an intermediate product in the DNA demethylation pathway, but may act as a functional epigenetic modification. We have conducted the largest study of site-specific 5hmC in placenta to date using parallel bisulphite and oxidative bisulphite modification with array-based assessment. Incorporating parallel RNA-sequencing data allowed us to assess associations between 5hmC and gene expression, using expression quantitative trait hydroxymethylation (eQTHM) analysis. We identified ~ 47,000 loci with consistently elevated (systematic) 5hmC proportions. Systematic 5hmC was significantly depleted (p < 0.0001) at CpG islands (CGI), and enriched (p < 0.0001) in 'open sea' regions (CpG >4 kb from CGI). 5hmC was most and least abundant at CpGs in enhancers and active transcription start sites (TSS), respectively (p < 0.05). We identified 499 significant (empirical-p <0.05) eQTHMs within 1 MB of the assayed gene. At most (75.4%) eQTHMs, the proportion of 5hmC was positively correlated with transcript abundance. eQTHMs were significantly enriched among enhancer CpGs and depleted among CpGs in active TSS (p < 0.05 for both). Finally, we identified 107 differentially hydroxymethylated regions (DHMRs, p < 0.05) across 100 genes. Our study provides insight into placental distribution of 5hmC, and sheds light on the functional capacity of this epigenetic modification in placenta.

Cell-Free DNA 5-Hydroxymethylcytosine Signatures for Lung Cancer Prognosis.

Accurate prognostic markers are essential for guiding effective lung cancer treatment strategies. The level of 5-hydroxymethylcytosine (5hmC) in tissue is independently associated with overall survival (OS) in lung cancer patients. We explored the prognostic value of cell-free DNA (cfDNA) 5hmC through genome-wide analysis of 5hmC in plasma samples from 97 lung cancer patients. In both training and validation sets, we discovered a cfDNA 5hmC signature significantly associated with OS in lung cancer patients. We built a 5hmC prognostic model and calculated the weighted predictive scores (wp-score) for each sample. Low wp-scores were significantly associated with longer OS compared to high wp-scores in the training [median 22.9 versus 8.2 months; p = 1.30 × 10-10; hazard ratio (HR) 0.04; 95% confidence interval (CI), 0.00-0.16] and validation (median 18.8 versus 5.2 months; p = 0.00059; HR 0.22; 95% CI: 0.09-0.57) sets. The 5hmC signature independently predicted prognosis and outperformed age, sex, smoking, and TNM stage for predicting lung cancer outcomes. Our findings reveal critical genes and signaling pathways with aberrant 5hmC levels, enhancing our understanding of lung cancer pathophysiology. The study underscores the potential of cfDNA 5hmC as a superior prognostic tool for guiding more personalized therapeutic strategies for lung cancer patients.

5-Hydroxymethylcytosine signals in serum are a predictor of chemoresistance in high-grade serous ovarian cancer.

OBJECTIVE: The genome-wide profiling of 5-hydroxymethylcytosines (5hmC) on circulating cell-free DNA (cfDNA) has revealed promising biomarkers for various diseases. The purpose of this study was to investigate 5hmC signals in serum cfDNA and identify novel predictive biomarkers for the development of chemoresistance in high-grade serous ovarian cancer (HGSOC). We hypothesized that 5hmC profiles in cfDNA reflect the development of chemoresistance and elucidate pathways that may drive chemoresistance in HGSOC. Moreover, we sought to identify predictors that would better stratify outcomes for women with intermediate-sensitive HGSOC. METHODS: Women diagnosed with HGSOC and known platinum sensitivity status were selected for this study. Nano-hmC-Seal was performed on cfDNA isolated from archived serum samples, and differential 5hmC features were identified using DESeq2 to establish a model predictive of chemoresistance. RESULTS: A multivariate model consisting of three features (preoperative CA-125, largest residual implant after surgery, 5hmC level of OSGEPL), stratified samples from intermediate sensitive, chemo-naive women diagnosed with HGSOC into chemotherapy-resistant- and sensitive-like strata with a significant difference in overall survival (OS). Independent analysis of The Cancer Genome Atlas data further confirmed that high OSGEPL1 expression is a favorable prognostic factor for HGSOC. CONCLUSIONS: We have developed a novel multivariate model based on clinico-pathologic data and a cfDNA-derived 5hmC modified gene, OSGEPL1, that predicted response to platinum-based chemotherapy in intermediate-sensitive HGSOC. Our multivariate model applies to chemo-naïve samples regardless if the patint was treated with adjuvant or neoadjuvant chemotherapy. These results merit further investigation of the predictive capability of our model in larger cohorts.

The Role of 5-Hydroxymethylcytosine as a Potential Epigenetic Biomarker in a Large Series of Thyroid Neoplasms.

Cytosine modifications at the 5-carbon position play a critical role in gene expression regulation and have been implicated in cancer development. 5-Hydroxymethylcytosine (5hmC), arising from 5-methylcytosine (5-mC) oxidation, has shown promise as a potential malignancy marker due to its depletion in various human cancers. However, its significance in thyroid tumors remains underexplored, primarily due to limited data. In our study, we evaluated 5hmC expression levels by immunohistochemistry in a cohort of 318 thyroid tumors. Our analysis revealed significant correlations between 5hmC staining extension scores and nodule size, vascular invasion, and oncocytic morphology. Nuclear 5hmC staining intensity demonstrated associations with focality, capsule status, extrathyroidal extension, vascular invasion, and oncocytic morphology. Follicular/oncocytic adenomas exhibited higher 5hmC expression than uncertain malignant potential (UMP) or noninvasive follicular thyroid neoplasms with papillary-like nuclear features (NIFTP), as well as malignant neoplasms, including papillary thyroid carcinomas (PTCs), oncocytic carcinomas (OCAs), follicular thyroid carcinomas (FTCs), and invasive encapsulated follicular variants of PTC (IEFV-PTC). TERT promoter mutation cases showed notably lower values for the 5hmC expression, while RAS (H, N, or K) mutations, particularly HRAS mutations, were associated with higher 5hmC expression. Additionally, we identified, for the first time, a significant link between 5hmC expression and oncocytic morphology. However, despite the merits of these discoveries, we acknowledge that 5hmC currently cannot segregate minimally invasive from widely invasive tumors, although 5hmC levels were lower in wi-FPTCs. Further research is needed to explore the potential clinical implications of 5hmC in thyroid tumors.

5-Hydroxymethylcytosine: the many faces of the sixth base of mammalian DNA.

Epigenetic phenomena play a central role in cell regulatory processes and are important factors for understanding complex human disease. One of the best understood epigenetic mechanisms is DNA methylation. In the mammalian genome, cytosines (C) in CpG dinucleotides were long known to undergo methylation at the 5-position of the pyrimidine ring (mC). Later it was found that mC can be oxidized to 5-hydroxymethylcytosine (hmC) or even further to 5-formylcytosine (fC) and to 5-carboxylcytosine (caC) by the action of 2-oxoglutarate-dependent dioxygenases of the TET family. These findings unveiled a long elusive mechanism of active DNA demethylation and bolstered a wave of studies in the area of epigenetic regulation in mammals. This review is dedicated to critical assessment of recent data on biochemical and chemical aspects of the formation and conversion of hmC in DNA, analytical techniques used for detection and mapping of this nucleobase in mammalian genomes as well as epigenetic roles of hmC in DNA replication, transcription, cell differentiation and human disease.

5-Hydroxymethylcytosine Profiling of Cell-Free DNA Identifies Bivalent Genes That Are Prognostic of Survival in High-Risk Neuroblastoma.

PURPOSE: Neuroblastoma is the most common extracranial solid tumor in childhood. We previously showed that circulating cell-free DNA (cfDNA) and tumor biopsy derived 5-hydroxymethylcytosime (5-hmC) profiles identified patients with neuroblastoma who experienced subsequent relapse. Here, we hypothesized that 5-hmC modifications selectively enriched in cfDNA compared with tumor biopsy samples would identify epigenetic changes associated with aggressive tumor behavior and identify novel biomarkers of outcome in patients with high-risk neuroblastoma. METHODS: 5-hmC profiles from cfDNA (n = 64) and tumor biopsies (n = 48) were compared. Two neuroblastoma cell lines underwent chromatin immunoprecipitation followed by sequencing (ChIP-Seq) for H3K27me3, H3K4me3, and H3K27ac; kethoxal-associated single-stranded DNA sequencing; hmC-Seal for 5-hmC; and RNA-sequencing (RNA-Seq). Genes enriched for both H3K27me3 and H3K4me3 in the included cell lines were defined as bivalent. Using bivalent genes defined in vitro, a bivalent signature was established in three publicly available cohorts of patients with neuroblastoma through gene set variation analysis. Differences between tumors with high or low bivalent signatures were assessed by the Kaplan-Meier method and Cox proportional hazards models. RESULTS: In cfDNA compared with tumor biopsy derived 5-hmC profiles, we found increased 5-hmC deposition on Polycomb Repressive Complex 2 target genes, a finding previously described in the context of bivalent genes. We identified 313 genes that bore bivalent chromatin marks, were enriched for mediators of neuronal differentiation, and were transcriptionally repressed across a panel of heterogeneous neuroblastoma cell lines. In three distinct clinical cohorts, low bivalent signature was significantly and independently associated with worse clinical outcome in patients with high-risk neuroblastoma. CONCLUSION: Low expression of bivalent genes is a biomarker of worse outcome in patients with high-risk neuroblastoma.

Targeting SIRT4/TET2 Signaling Alleviates Human Keratinocyte Senescence by Reducing 5-hydroxymethylcytosine Loss.

Skin aging is characterized by wrinkle formation and increased frailty and laxity, leading to the risk of age-related skin diseases. Keratinocyte is an important component of the epidermis in skin structure, and keratinocyte senescence has been identified as a pivotal factor in skin aging development. Because epigenetic pathways play a vital role in the regulation of skin aging, we evaluated human skin samples for DNA hydroxymethylation (5-hydroxymethylcytosine; 5-hmC) and SIRT4 expressions. Results found that both 5-hmC and SIRT4 showed a significant decrease in aged human skin samples. To test the results in vitro, human keratinocytes were cultured in H2O2, which modulates skin aging in vivo. However, H2O2-induced keratinocytes showed senescence-associated protein expression and significant downregulation of 5-hmC and SIRT4 expressions. Moreover, 5-hmC-converting enzymes ten eleven translocation 2 (TET2) showed a decrease and enhanced TET2 acetylation level in H2O2-induced keratinocytes. However, the overexpression of SIRT4 in keratinocytes alleviates the senescence phenotype, such as senescence-associated protein expression, decreases the TET2 acetylation, but increases TET2 and 5-hmC expressions. Our results provide a novel relevant mechanism whereby the epigenetic regulation of keratinocytes in skin aging may be correlated with SIRT4 expression and TET2 acetylation in 5-hmC alteration. Our study may provide a potential strategy for antiskin aging, which targets the SIRT4/TET2 axis involving epigenetic modification in keratinocyte senescence.

Locus-Specific Enrichment Analysis of 5-Hydroxymethylcytosine Reveals Novel Genes Associated with Breast Carcinogenesis.

An imbalance in DNA methylation is a hallmark epigenetic alteration in cancer. The conversion of 5-methylcytosine (5-mC) to 5-hydroxymethyl cytosine (5-hmC), which causes the imbalance, results in aberrant gene expression. The precise functional role of 5-hydroxymethylcytosine in breast cancer remains elusive. In this study, we describe the landscape of 5-mC and 5-hmC and their association with breast cancer development. We found a distinguishable global loss of 5-hmC in the localized and invasive types of breast cancer that strongly correlate with TET expression. Genome-wide analysis revealed a unique 5-mC and 5-hmC signature in breast cancer. The differentially methylated regions (DMRs) were primarily concentrated in the proximal regulatory regions such as the promoters and UTRs, while the differentially hydroxymethylated regions (DhMRs) were densely packed in the distal regulatory regions, such as the intergenic regions (>-5 kb from TSSs). Our results indicate 4809 DMRs and 4841 DhMRs associated with breast cancer. Validation of nine 5-hmC enriched loci in a distinct set of breast cancer and normal samples positively correlated with their corresponding gene expression. The novel 5-hmC candidates such as TXNL1, and CNIH3 implicate a pro-oncogenic role in breast cancer. Overall, these results provide new insights into the loci-specific accumulation of 5-mC and 5-hmC, which are aberrantly methylated and demethylated in breast cancer.

Regional gain and global loss of 5-hydroxymethylcytosine coexist in genitourinary cancers and regulate different oncogenic pathways.

BACKGROUND: DNA 5-hydroxymethylcytosine (5hmC) is produced by dynamic 5mC oxidation process contributing to tissue specification, and loss of 5hmC has been reported in multiple cancers including genitourinary cancers. However, 5hmC is also cell-type specific, and its variability may exist between differentiated tumor cells and cancer stem cells. Thus, cancer-associated changes in 5hmC may be contributed by distinct sets of tumor cells within the tumor tissues. RESULTS: Here, we applied a sensitive immunoprecipitation-based method (hMeDIP-seq) to analyze 5hmC changes during genitourinary carcinogenesis (including prostate, urothelial and kidney). We confirmed the tissue-specific distribution of 5hmC in genitourinary tissues and identified regional gain and global loss of 5hmC coexisting in genitourinary cancers. The genes with gain of 5hmC during tumorigenesis were functionally enriched in regulating stemness and hypoxia, whereas were associated with poor clinical prognosis irrespective of their differences in tumor type. We identified that gain of 5hmC occurred in soft fibrin gel-induced 3D tumor spheres with a tumor-repopulating phenotype in two prostate cancer cell lines, 22RV1 and PC3, compared with conventional two-dimensional (2D) rigid dishes. Then, we defined a malignant signature derived from the differentially hydroxymethylated regions affected genes of cancer stem-like cells, which could predict a worse clinical outcome and identified phenotypically malignant populations of cells from prostate cancer tumors. Notably, an oxidation-resistant vitamin C derivative, ascorbyl phosphate magnesium, restored 5hmC and killed the cancer stem cell-like cells leading to apoptosis in prostate cancer cell lines. CONCLUSIONS: Collectively, our study dissects the regional gain of 5hmC in maintaining cancer stem-like cells and related to poor prognosis, which provides proof of concept for an epigenetic differentiation therapy with vitamin C by 5hmC reprogramming.

Genome-wide profiling of 5-hydroxymethylcytosines in circulating cell-free DNA reveals population-specific pathways in the development of multiple myeloma.

Multiple myeloma (MM) and its precursors monoclonal gammopathy of undetermined significance (MGUS) and smoldering myeloma (SMM) are 2-3 times more common in African Americans (AA) than European Americans (EA). Although epigenetic changes are well recognized in the context of myeloma cell biology, the contribution of 5-hydroxymethylcytosines (5hmC) to racial disparities in MM is unknown. Using the 5hmC-Seal and next-generation sequencing, we profiled genome-wide 5hmC in circulating cell-free DNA (cfDNA) from 342 newly diagnosed patients with MM (n = 294), SMM (n = 18), and MGUS (n = 30). We compared differential 5hmC modifications between MM and its precursors among 227 EA and 115 AA patients. The captured 5hmC modifications in cfDNA were found to be enriched in B-cell and T-cell-derived histone modifications marking enhancers. Of the top 500 gene bodies with differential 5hmC levels between MM and SMM/MGUS, the majority (94.8%) were distinct between EA and AA and enriched with population-specific pathways, including amino acid metabolism in AA and mainly cancer-related signaling pathways in EA. These findings improved our understanding of the epigenetic contribution to racial disparities in MM and suggest epigenetic pathways that could be exploited as novel preventive strategies in high-risk populations.

iRNA5hmC-HOC: High-order correlation information for identifying RNA 5-hydroxymethylcytosine modification.

RNA 5-hydroxymethylcytosine (5 hmC) is an important RNA modification, which plays vital role in several biological processes. Currently, it is a hot topic to identify 5 hmC sites due to its benefit in understanding its biological functions. Therefore, in this study, we developed a predictor called iRNA5 hmC-HOC, which is based on a high-order correlation information method to identify 5 hmC sites. To build the model, 22 different classes of dinucleotide physicochemical (PC) properties were employed to represent RNA sequences, and the least absolute shrinkage and selection operator (LASSO) algorithm was adopted to select the most discriminative features. In the jackknife test, the proposed method achieved 89.80% classification accuracy based on support vector machine (SVM). As compared with the state-of-the-art predictors, our proposed method has significant improvement on the classification performance. It indicates that the proposed method might be a promising tool in identifying RNA 5 hmC modification sites. The dataset and source codes are available at https://figshare.com/articles/online_resource/iRNA5hmC-HOC/15177450.

R5hmCFDV: computational identification of RNA 5-hydroxymethylcytosine based on deep feature fusion and deep voting.

RNA 5-hydroxymethylcytosine (5hmC) is a kind of RNA modification, which is related to the life activities of many organisms. Studying its distribution is very important to reveal its biological function. Previously, high-throughput sequencing was used to identify 5hmC, but it is expensive and inefficient. Therefore, machine learning is used to identify 5hmC sites. Here, we design a model called R5hmCFDV, which is mainly divided into feature representation, feature fusion and classification. (i) Pseudo dinucleotide composition, dinucleotide binary profile and frequency, natural vector and physicochemical property are used to extract features from four aspects: nucleotide composition, coding, natural language and physical and chemical properties. (ii) To strengthen the relevance of features, we construct a novel feature fusion method. Firstly, the attention mechanism is employed to process four single features, stitch them together and feed them to the convolution layer. After that, the output data are processed by BiGRU and BiLSTM, respectively. Finally, the features of these two parts are fused by the multiply function. (iii) We design the deep voting algorithm for classification by imitating the soft voting mechanism in the Python package. The base classifiers contain deep neural network (DNN), convolutional neural network (CNN) and improved gated recurrent unit (GRU). And then using the principle of soft voting, the corresponding weights are assigned to the predicted probabilities of the three classifiers. The predicted probability values are multiplied by the corresponding weights and then summed to obtain the final prediction results. We use 10-fold cross-validation to evaluate the model, and the evaluation indicators are significantly improved. The prediction accuracy of the two datasets is as high as 95.41% and 93.50%, respectively. It demonstrates the stronger competitiveness and generalization performance of our model. In addition, all datasets and source codes can be found at https://github.com/HongyanShi026/R5hmCFDV.

5-Hydroxymethylcytosine profiles in plasma cell-free DNA reflect molecular characteristics of diabetic kidney disease.

Background: Diabetic kidney disease (DKD), one of the main complications of diabetes mellitus (DM), has become a frequent cause of end-stage renal disease. A clinically convenient, non-invasive approach for monitoring the development of DKD would benefit the overall life quality of patients with DM and contribute to lower medical burdens through promoting preventive interventions. Methods: We utilized 5hmC-Seal to profile genome-wide 5-hydroxymethylcytosines in plasma cell-free DNA (cfDNA). Candidate genes were identified by intersecting the differentially hydroxymethylated genes and differentially expressed genes from the GSE30528 and GSE30529. Then, a protein interaction network was constructed for the candidate genes, and the hub genes were identified by the MCODE and cytoHubba algorithm. The correlation analysis between the hydroxymethylation level of the hub genes and estimated glomerular filtration rate (eGFR) was carried out. Finally, we demonstrated differences in expression levels of the protein was verified by constructing a mouse model of DKD. In addition, we constructed a network of interactions between drugs and hub genes using the Comparative Toxicogenomics Database. Results: This study found that there were significant differences in the overall distribution of 5hmC in plasma of patients with DKD, and an alteration of hydroxymethylation levels in genomic regions involved in inflammatory pathways which participate in the immune response. The final 5 hub genes, including (CTNNB1, MYD88, CD28, VCAM1, CD44) were confirmed. Further analysis indicated that this 5-gene signature showed a good capacity to distinguish between DKD and DM, and was found that protein levels were increased in renal tissue of DKD mice. Correlation analysis indicated that the hydroxymethylation level of 5 hub genes were nagatively correlated with eGFR. Toxicogenomics analysis showed that a variety of drugs for the treatment of DKD can reduce the expression levels of 4 hub genes (CD44, MYD88, VCAM1, CTNNB1). Conclusions: The 5hmC-Seal assay was successfully applied to the plasma cfDNA samples from a cohort of DM patients with or without DKD. Altered 5hmC signatures indicate that 5hmC-Seal has the potential to be a non-invasive epigenetic tool for monitoring the development of DKD and it provides new insight for the future molecularly targeted anti-inflammation therapeutic strategies of DKD.

Cell-free DNA 5-hydroxymethylcytosine is an emerging marker of acute myeloid leukemia.

Aberrant changes in 5-hydroxymethylcytosine (5hmC) are a unique epigenetic feature in many cancers including acute myeloid leukemia (AML). However, genome-wide analysis of 5hmC in plasma cell-free DNA (cfDNA) remains unexploited in AML patients. We used a highly sensitive and robust nano-5hmC-Seal technology and profiled genome-wide 5hmC distribution in 239 plasma cfDNA samples from 103 AML patients and 81 non-cancer controls. We developed a 5hmC diagnostic model that precisely differentiates AML patients from controls with high sensitivity and specificity. We also developed a 5hmC prognostic model that accurately predicts prognosis in AML patients. High weighted prognostic scores (wp-scores) in AML patients were significantly associated with adverse overall survival (OS) in both training (P = 3.31e-05) and validation (P = 0.000464) sets. The wp-score was also significantly associated with genetic risk stratification and displayed dynamic changes with varied disease burden. Moreover, we found that high wp-scores in a single gene, BMS1 and GEMIN5 predicted OS in AML patients in both the training set (P = 0.023 and 0.031, respectively) and validation set (P = 9.66e-05 and 0.011, respectively). Lastly, our study demonstrated the genome-wide landscape of DNA hydroxymethylation in AML and revealed critical genes and pathways related to AML diagnosis and prognosis. Our data reveal plasma cfDNA 5hmC signatures as sensitive and accurate markers for AML diagnosis and prognosis. Plasma cfDNA 5hmC analysis will be an effective and minimally invasive tool for AML management.